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Abstract 

We perform a novel type of analysis of diffractive deep-inelastic scattering data, in 



turbative QCD expressions. In particular, we treat individually the components of the 
Pomeron of different size. We are able to describe simultaneously both the recent ZEUS 
"q^' and HI diffractive data. In addition to the usual two-gluon model for the perturbative 

o ■ Pomeron, we allow for the possibility that it may be made from two sea quarks. 

D ; 

"^ ' A notable feature of deep-inelastic scattering is the existence of diffractive events, 7*p -^ Xp, 

in which the slightly deflected proton and the cluster X of outgoing hadrons are well-separated 

rS [ in rapidity. The large rapidity gap is believed to be associated with Pomeron exchange. The 
a ' diffractive events make up an appreciable fraction of all (inclusive) deep-inelastic events, 7*p — *■ 
X. We will refer to the diffractive and inclusive processes as DDIS and DIS respectively. 

Here we perform a perturbative QCD analysis of the new high precision DDIS data, recently 
obtained by the ZEUS 012] and HI [3] Collaborations at HERA. The analysis is novel in that it 
treats individually the components of the Pomeron of different transverse size. The description 
of the DDIS data is based on a purely perturbative QCD framework. We take input forms of the 
parton distributions of the Pomeron given by the calculation of the lowest-order QCD diagrams 
for 7*p -^ Xp Q. In previous analyses, the Pomeron was treated as a hadron-like object of 
more or less flxed size. However, the microscopic structure of the Pomeron is different to that 
of a hadron. In perturbative QCD, it is known that the Pomeron singularity is not an isolated 
pole, but a branch cut, in the complex angular momentum plane [Hj. The pole singularity 
corresponds to a single particle, whereas a branch cut may be regarded as a continuum series 
of poles. That is, the Pomeron wave function consists of a continuous number of components. 
Each component i has its own size, l//ij. The QCD DGLAP evolution of a component should 
start from its own scale /Xj, provided that //« is large enough for the perturbative evolution to be 
valid. Therefore, the expression for the diffractive structure function F.^ contains an integral 
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Figure 1: (a) Cut diagram illustrating the main idea contained in (^. For each component of the 
perturbative Pomeron of size l//i, represented by two t-channel gluons in a colour singlet, the 
Pomeron structure function -Ff (/3, Q"^; /i^) is evaluated from the quark singlet, [3TF{(3, Q^\ /i^), 
and gluon, (3(f^{(5, Q^\ /z^), distributions of the Pomeron. The perturbative Pomeron flux factor 
fip{xip; /i^) is given in terms of the gluon distribution of the proton, xp g{xip, iJ?). (b) Later, we 
also include contributions from diagrams in which the Pomeron is represented by sea quark- 
antiquark exchange (plus interference with the two-gluon Pomeron). 



over the Pomeron size, or rather over the scale /i. So to obtain F^ we evolve the input parton 
distributions of each component of the Pomeron from their own starting scale /x up to the final 
scale Q. The extra integral over /x reflects the fact that the partonic structure of the Pomeron 
is more complicated than that of a normal hadron. 

Recall that in the usual analyses (for example, by HI P|) Regge factorisation P is as- 
sumed, such that the diffractive structure function Fg (x^p,/?, Q^, t)^ is written as a product 
of the Pomeron flux factor fip{xjp,t) and the structure function F^{[3,Q'^) which describes 
the interaction of the Pomeron with the virtual photon probe. The input Pomeron parton 
distributions are taken to be arbitrary polynomials. The Pomeron flux factor is taken from 
Regge phenomenology with some effective Pomeron intercept aip(O). However, the value of 
ttip(O) = 1.17 [2] needed to fit DDIS data is significantly higher than the value of 1.08 obtained 
from soft hadron data [7j. Instead, in the present analysis, we have /x-factorisation, such that 
the t-integrated observable is 



F?^^^\xjp,P,Q' 



dfi^ fjp{xjp;fi^)F^{P,Q^-fi' 



(1) 



see Fig. ^ Here, the subscript P on Fg p is to indicate that this is the perturbative contribution 

with yU > Qo ~ 1 GeV. We will introduce the full form which we take for F2 in a moment. 

The Pomeron structure function, Ff (/?, Q^; /i^), is obtained by NLO DGLAP evolution up 

to Q^ from input Pomeron parton distributions parameterised at a starting scale /x^. The 

^Here, xip is the fraction of the proton's momentum transferred through the rapidity gap by the Pomeron, 
/3 = Xb/xip is the fraction of the Pomeron's momentum carried by the struck quark, Xb is the Bjorken x 
variable, Q^ is the photon virtuality, and t is the squared 4-momentum transfer. 





Figure 2: (a) Quark dipole and (b) effective gluon dipole interacting with the proton via a 
perturbative Pomeron composed of two t-channel gluons. Here, l± is a space-hke 4-vector such 
that /5_ = ~^t) and fg{x]p,lf, fi'^) is the unintegrated gluon distribution of the proton. 



P dependence of these input distributions is obtained from lowest-order perturbative QCD 
diagrams. The perturbative Pomeron is represented by two t-channel gluons in a colour singlet; 
we will indicate this using the notation IP = G. The relevant lowest-order diagrams are 
shown in Fig. |21 We regard the virtual photon as dissociating into either a quark-ant iquark 
dipole (Fig. I2ta)) or an effective gluon dipole, made up of a gluon and a compact qq pair 
(Fig. Efb)). In each case the four permutations of the couplings of the two t-channel gluons to 
the two components of the dipole are implied. Evaluating these diagrams using light-cone wave 
functions of the photon, in the limit of strongly-ordered transverse momenta {k '^ kt <^ Q), 
leads to the input quark singlet and gluon distributions of the form jlj 

/?S^=^(/3,/x2;/x2) = c,/g/3^(1-/3), (2) 

;i + 2/3')' (1 - /3')^ (3) 



/3r=^(/3',/z2;/i^) 



-alG 



where (3' is the Pomeron's momentum fraction carried by the gluon with momentum k in 
Fig. I21^b). The parameters Cq/c and Cg/c implicitly include all the numerical factors arising 
from the lowest-order calculations. We will let these normalisations go free in fits to the DDIS 
data to account for higher-order QCD corrections (effective i^-factors). The quark singlet 
distribution is TF ^u + d + s + u + d + s, with u = d = s = u = d = s,so that the non-singlet 
distributions are all zero. The contributions of the charm and bottom quarks to F![{(3, Q^\ /i^) 
are calculated in the heavy quark fixed-flavour number scheme. 

Moreover, for the perturbative contribution, the Pomeron flux factor jjp[xjp] /i^) is given in 
terms of the integrated gluon distribution of the proton. Consider the lowest-order diagram. 
Fig. |21^a), in which the virtual photon dissociates into a forward-going quark-antiquark pair, 
that is, 7* — i> gg as t ^ 0J |H] . After summation over the four amplitudes, the integral of the 
unintegrated gluon distribution /^(xjp, /^,/i^) over the gluon loop transverse momentum l\ up 



to the quark virtuality /i^ = k^lil — [3) gives the integrated gluon distribution xjp g{xjp, fi'^). 
Thus the Pomeron flux factor is 



flP=G{.Xip] i/) = — 
Xjp 



2^ 

(4) 



«5(/i ; . 2 
; — Xipg{xip,iJ 



Here, a factor 1/Bd from the t-integration, where B^, ~ 6 GeV^^ is the diffractive slope pa- 
rameter, has been absorbed into the parameters Cg/c and Cg/c of (JD) and Q. Strictly speaking, 
the Pomeron flux factor (j^ should be written in terms of the skewed gluon distribution. At 
small Xjp this gives rise to an overall constant factor |9,, i?^, which again we absorb into the 
parameters Cg/c and Cg/c- 

In addition to the leading-twist contribution arising from Fig. |2ta) with a transversely 
polarised photon, there is an analogous twist-four contribution to F2 arising from Fig. Efa) 
with a longitudinally polarised photon, 

F^^^\xjp,(3,Q') = n^ df,' ^ M^ip-.f^')] FfiP). (5) 

The (3 dependence is again obtained from lowest-order perturbative QCD calculations |3]: 

i^f=''(/5) = CL/G/5^(2/3-l)^ (6) 

where, as before, cl/g is taken to be a free parameter. The twist-four nature of this longitudinal 
contribution is evident from the fi^/Q"^ factor in (0). 

We also include a non-perturbative (NP) Pomeron contribution (from scales fi < Qq) and a 
secondary Reggeon (IR) contribution to Fg (xjp,/?, Q^), so that 

pD(3) _ pD(3) pD{3) pD(3) pZ){3) /7^ 

with 

f3?(x^, /?, Q2) = /^=Np(xip) Ff =^P(/?, Q^ g2)^ (8) 

F^g^ixjp, A g2) = cjR Uxjp) Ff (/5, Q2), (9) 

where qr is taken to be a free parameter. Here, the non-perturbative Pomeron and Reggeon 
flux factors are^ 

l-2ai(0) f^ _ Bacut -20itcut\ 

/'t'lnin pi'i ^ fP \ IP I 

with i = IP and Ji? respectively, and aj(t) = q;j(0) + a^ t. The integration limits are taken to be 
^cut = —1 GeV^ and tmin ~ GeV^. For the non-perturbative Pomeron, we flx ajp{0) = 1.08 
[7j, a'jp = 0.26 GeV~^, and Bjp = 4.6GeV~^ [TU], whereas for the secondary Reggeon we take 
affi(O) = 0.50 [nj, a^ = 0.90 GeV"^ [12], and Bm = 2.0 GeV'^ ^. Apart from aip{0), 

^The couplings of the Pomeron or secondary Reggeon to the proton are absorbed into the parameters Cq/np 
(a = q,g) and qr. 



these are the same values used in the prehminary HI analysis [3]. The secondary Reggeon 
structure function, F^(/?, Q^), is calculated at NLO from the GRV pionic parton distribu- 
tions ^^. For the non-perturbative Pomeron, the input quark singlet and gluon distributions, 
PTF"^^^ {P , Ql; Ql) and P' (/^^^^ {P' , Qq', Qq) , are taken to have the same P dependence as for 
the two-quark Pomeron introduced later on (see ()14|) and (jlSp ). with different normalisations 
Cq/NP and Cg/Np. (Taking the same P dependence as for the two-gluon Pomeron, and Q, 
gives a much worse description of the data.) 

We fit to the preliminary ZEUS 0121 and HI |3] DDIS data using ((Tj), and varying the free 
parameters until an optimum description of the data is obtained. We impose a cut Mx > 2 
GeV on the fitted data to exclude large contributions from vector meson production and other 
higher- twist effects, and a cut y < 0.45 so that we can assume that the measured reduced 
diffractive cross section cTr- is approximately equal to F2 . The statistical and systematic 
experimental errors are added in quadrature. We use the qcdnum program ^3] to perform 
the NLO DGLAP evolution and the minuit program |16j to find the optimal parameters. The 
values of as{M'^) and the charm and bottom quark masses are taken to be the same as in the 
MRST2001 NLO parton set ^7). Two sets of preliminary ZEUS data are fitted: those obtained 
using the leading proton spectrometer (LPS) [T], and those obtained using the so-called Mx 
method [2] which is based on the fact that diffractive and non-diffractive events have very 
different InMjr distributions. For the latter data set, in addition to elastic proton scattering, 
proton dissociation up to mass My = 2.3 GeV is included. Clearly the cross section will be 
larger in this case, so we allow for the overall normalisation of these data by multiplying ((Tj) 
by a factor Nz- An analogous normalisation, Nh, is applied for the preliminary HI data [3], 
where diffractive events are selected on the basis of a large rapidity gap, and where proton 
dissociation up to mass My = 1.6 GeV is included. The ZEUS Mx data i2j do not include 
secondary Reggeon contributions, therefore we omit the fourth term of (I7|) when fitting to 
these data. We fit to each data set separately, and then we perform fits to the three data sets 
combined. 

For our first study, we parameterise the perturbative Pomeron fiux factor (0)) using a sim- 
plified form for the gluon distribution of the proton, 

xjpg{xjp,f/) =Xjp^, (11) 

where A is independent of /i^ and is determined by the fit to data.^ The normalisation of (jll|) 
has been absorbed into the free parameters Cq/c, Cg/c, and cl/g- 

Varying the Qo parameter, we find that the best fit to the combined ZEUS and HI data 
sets is obtained with Ql = 0.8 GeV^, which gives a x^/d.o.f. = 1.05 with c^/np going to zero. 
Later on, we will use the MRST2001 NLO ^7j parton distributions of the proton instead of 
the simplified form ()lip. where the minimum possible scale is 1 GeV. Using the form ()11|) with 

^Strictly speaking, A should depend on In^^. We investigated this effect by talcing A(^^) = 0.08 + 
CAln(/xV(0.45GeV^)) with Qo = 1 GeV and c^/np = 0. The combined fit to ZEUS and HI DDIS data 
gave a x^/d.o.f. = 1.12 with c\ ~ 0.054 ± 0.006. This is consistent with the value found by HI in a fit to 
inclusive F2 data [iBj of ca = 0.0481 ± 0.0013(stat.) ± 0.0037(syst.). Since the x^/d.o.f. was not improved 
compared to the corresponding fit which took A to be independent of fi^ (x^/d.o.f. = 1.07), we used the form 
ifTTI for simplicity. 



Data sets fitted 


ZEUS LPS 


ZEUS Mx 


HI 


ZEUS ± HI 


Number of points 


69 


121 


214 


404 


xVd.o.f. 


0.67 


0.78 


1.08 


1.08 


c,iG (GeV2) 


0.71 ±0.39 


0.48 ±0.12 


2.2 ±0.4 


1.13 ±0.15 


c,,G (GeV2) 


0.11 ±0.05 


0.10 ±0.02 


0.26 ±0.05 


0.17±0.02 


CLIG (GeV2) 





0.20 ±0.08 


0.54 ±0.17 


0.36 ±0.08 


c,/NP (GeV-2) 


0.87 ±0.13 


1.22 ±0.04 


0.91 ±0.05 


1.09 ±0.05 


cjR (GeV-2) 


6.7 ±0.8 


— 


7.5 ±2.0 


6.2 ±0.6 


A 


0.23 ±0.04 


0.21 ±0.02 


0.13 ±0.01 


0.17±0.01 


iVz 




1.56 (fixed) 


— 


1.56 ±0.06 


iV^ 




— 


1.26 (fixed) 


1.26 ±0.05 


i?(6.5GeV2),i?(90GeV2) 


0.60,0.60 


0.56,0.57 


0.54,0.55 


0.55,0.56 



Table 1: Tlie values of the free parameters obtained in the fits to preliminary ZEUS ^1121 and 
HI 0-^2 data with a gluon distribution of the proton proportional to x'^ ()11|) . The last 
row RlQ"^), defined in (J2H), gives the fraction of the Pomeron's (plus Reggeon's) momentum 
carried by gluons at xp = 0.003. 

Ql = 1 GeV^ gives only a slightly worse x^/d.o.f. = 1.07. Furthermore, fixing c^/np = makes 
little difference to the quality of the fit. Therefore, in all fits presented here, we take Qq = 1 
GeV^ and fix c^/np = 0. 

We find that each data set can be well described by this simple, perturbatively-motivated, 
approach. However, different values of A and the other parameters are obtained from the ZEUS 
and HI data, as can be seen from Table ^ In particular, the HI data seem to have a flatter 
X]p dependence than the ZEUS data. This should be regarded as some inconsistency between 
the data sets, but not as a contradiction, since it is possible to obtain an adequate description 
of the combined data sets, as shown in Fig. |S1 and by the results in the last column of Table [T] 

These fits to the DDIS data imply that the growth of F2 with decreasing xjp comes from 
a gluon distribution which grows as x]^^ with A ~ 0.17. On the other hand, at low scales 
/x ~ Qo ~ 1 GeV, which are dominant due to the l//i^ factor in the Pomeron flux factor (JU, 
the gluon distribution of the proton obtained from global analyses of DIS and related data is 
valence-like, or even negative, at small x, while the sea quark distribution grows as a negative 
power of X with decreasing x. Therefore, in order to describe the DDIS data we are forced to 
introduce a Pomeron comprised of two t-channel sea quarks, illustrated in Fig. ^b). In analogy 
to the flux factor (^ for the two-gluon component, we therefore define a perturbative Pomeron 



1997 ZEUS LPS data (prel.) 

MRW 2004 NLO QCD fit 



Perturbative contrib. 
Non-perturbative contrib. 
Twist-4 contrib. 
Reggeon contrib. q= (Cev' 



1998/99 ZEUS data (prel.) 

MRW 2004 NLO QCD fit 



Perturbative contrib. 
Non-perturbative contrib. 
Twist-4 contrib. 



£2- 004 
Q CM 

u. 

0.02 

X 



P = 0.007 p = 0.030 p = 0.130 





q' (Gev') 



39 0.05 



1 1 1 1 1 1 1 1 1 1 1 1 



10 10 10 10 10 10 10 10 10 10 10 10 



•^IP 



MP 



1997 H1 data (prel.) 

MRW 2004 NLO QCD frt 



Perturbative contrib. 
Non-perturbative contrib. 



Tw/ist-4 contrib. 
1 contrib. 




10" W W icr" 10"^" W io" W 10^ -tr W io' to" To^ trir ^o" io 

X|p 



Figure 3: Fit to combined preliminary ZEUS ^1 |2j and HI \Z^ Fj data with a gluon distri- 
bution of the proton proportional to Xjp (fTTj) . The curves show the four contributions to the 
total, as defined in ((7j). Only data points included in the fit are plotted. 
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Figure 4: (a) Quark dipole and (b) effective gluon dipole interacting with the proton via 
a perturbative Pomeron composed of two t-channel sea quarks. Here, /^(xjp, /^,/x^) is the 
unintegrated sea quark distribution of the proton. 



flux factor for this two-quark component {IP = S) and for the interference term {IP = GS): 

(12) 



fip=s{xip;n ) 
fip=Gs{xip; /i^) 



1 

Xjp 

1 

Xjp 



-^ Xjp S{Xjp, /i^) 



as{fJ^ 



2 xjpg{xip, fi^) Xjp S{xip, fx^ 



(13) 



where S{xjp,fi'^) = 2[u{xjp, fi^) + d{xjp,^'^) + s{xjp,^'^)] is the integrated sea quark distribution 
of the proton. Just as for the starting distributions of the two-gluon Pomeron, (j2)) and Q, 
we calculate the jS dependence of the input Pomeron parton distributions for the two-quark 
Pomeron (and the interference contribution) from the diagrams in Fig. 0] using lowest-order 
perturbative QCD. Unlike for the two-gluon Pomeron, there are only two permutations of the 
couplings of the two t-channel sea quarks to the two components of the dipole. We find"^ 

pj:^='{P,fx';fx') = c,/sP{l-P), 



(14) 



-9/S' 



/3' 



7n2 



(3'g^='{(3',^^'■,^^') 
/?E^=^^(A/x2;/i2) = c,/G5/5'(l-/5), 



(15) 
(16) 

/3'/'=^^(/?', /x^; /x^) = c,/as (1 + 2/?') (1 - [3')'. (17) 

The f5 dependent factors of the twist-four contribution arising from Fig. El^a) with a longitu- 
dinally polarised photon (and the interference contribution), analogous to ©, are found to 
be 



7iP=5 



-,JP=GS 



iP) 

{(3) 



CL/GS (3^ (2/3 



(18) 
(19) 



*Note that (iYi 



'IP=S 



{P , jj? ; li^) has the same (3 dependence as the F^ term in the BEKW model [T^. 
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Data sets fitted 


ZEUS LPS 


ZEUS Mx 


HI 


ZEUS ± HI 


Number of points 


69 


121 


214 


404 


xVd.o.f. 


0.79 


0.96 


0.71 


1.14 


c,iG (GeV^) 


0.001 ±0.053 


0.018 ±0.023 


0.36 ±0.06 


0.18 ±0.04 


c,iG (GeV^) 








0.37 ±0.02 





CLIG (GeV^) 


0.21 ± 1.48 


0.050 ±0.033 


0.14 ±0.03 


0.064 ±0.024 


c,is (GeV2) 


0.97 ±0.40 


0.49 ±0.10 


1.06 ±0.13 


0.58 ±0.07 


c,is (GeV2) 


1.23 ±0.21 


1.23 ±0.07 





1.31 ±0.07 


CLis (GeV^) 


0.41 ±0.28 


0.21 ±0.09 





0.11 ±0.05 


c,/NP (GeV-2) 


0.79 ±0.22 


1.16 ±0.08 


0.09 ±0.11 


0.92 ±0.07 


CK (GeV-2) 


6.6 ±0.7 




8.4 ±1.8 


6.4 ±0.5 


^z 


— 


1.54 (fixed) 


— 


1.54 ±0.06 


A^H 


— 




1.24 (fixed) 


1.24 ±0.04 


i?(6.5GeV^),i?(90GeV2) 


0.57, 0.58 


0.57, 0.59 


0.60, 0.66 


0.57, 0.57 



Table 2: The values of the free parameters obtained in the fits to ZEUS ^1121 and HI [Sj Fj 
data with MRST2001 NLO [TTj gluon and sea quark distributions of the proton. The last row 
R{Q^), defined in ()21|). gives the fraction of the Pomeron's (plus Reggeon's) momentum carried 
by gluons at xjp = 0.003. 

The normalisation of the interference terms between the two-gluon and the two-quark Pomerons 
is fixed by Ci/cs = J^^ijG Ci/Si where i = q,g, L; that is, the i^-factor is fixed for the amplitude 
rather than for the cross section. 

The results of fits with this extended model, using the MRST2001 NLO fTTJ gluon and 
sea quark distributions of the proton, are shown in Table |21 We set xjp g{xip, fi'^) = if it is 
negative. Again, good fits are obtained whether fitting ZEUS and HI data separately or all 
together. However, the fit with only HI data is dramatically different from the other three fits 
in Table with a much larger two-gluon Pomeron contribution compared to the other three, 
which are dominated by the two-quark Pomeron. This difference can be traced to the flatter 
xip dependence of the HI data compared to the ZEUS data (see Table Q). Note that some 
parameters in Table |21 are consistent with zero, indicating some redundancy in this extended 
model. 

From these flts to Fg , we can extract diffractive parton distributions, a^{xjp, P,Q'^) = 
PT,^{xjp, P,Q'^) or j3g^{xjp, j3,Q^), from the three leading-twist contributions to ((Tj): 



TD /-"* o r^ a \ '^ (J.i\ 



IP=G,S,GS \-"^o 

+ fip=Np{xip) a^=^P(/?, Q^ Ql) + cjR Uxjp) a^{(3, Q^). (20) 



The diffractive parton distributions calculated using (pUj) are plotted for the fits to the combined 
ZEUS and HI data sets of Tables O ('A') and El ('MRST') in Fig. El^a) for xp = 0.003 and 
Q^ = 6.5,90 GeV^. We also show the Pomeron parton distributions from the preliminary HI 
analysis [3^ multiplied by fip{xip) (given by (fTUIl with aip{0) = 1.173) and normalised to the 
ZEUS LPS data by dividing by a factor 1.26 (from Table Q). The eight different fits of Tables 
n and 121 are found to give similar diffractive parton distributions, especially at the higher Q^ 
value, with the possible exception of the 'MRST' fit to only HI data. In Fig. Efb) we show the 
breakdown of the five separate components of (j^UI) for the 'MRST' fit to the combined data 
sets. Note the large contribution from the two-quark component of the Pomeron. 

From Fig. the diffractive quark singlet distribution obtained by HI 3J has a slightly 
steeper Q^ dependence than the fits presented here, and hence HI obtain a larger diffractive 
gluon distribution. In addition, the smaller value of ^^(Ml) used by HI also enlarges their 
gluon density.^ In our analysis, all the input Pomeron parton distributions vanish as either 
(1 — /5) or (1 — /5)^ as /5 — >■ 1. As /? — >■ 1, the only non-zero contribution to Fg comes from 
a twist-four component arising from longitudinally polarised photons. This contribution was 
not included in the HI analysis Ej, and hence rather large diffractive parton distributions were 
obtained by HI for j3 close to 1, with an unphysical 'bump' in the diffractive gluon distribution 
(see Fig. EI). 

We define the fraction of the Pomeron's (plus Reggeon's) momentum carried by gluons at 
xjp = 0.003 as 

^.Q2^ ^ J^,,dppg^{xjp = 0.003, P,Q') 

J^^^dp [PED{xjp = 0.003, P,Q'')+f3gD{xjp = 0.003, P,Q'')]' 

which is given for Q^ values of 6.5 and 90 GeV^ in the last rows of Tables Q] and |21 The 
gluon momentum fraction, R{Q'^), is consistently 55-60% and is almost independent of Q^. 
Taking the same as'(Ml) as in the preliminary HI analysis would increase this value to ~ 65%, 
compared to the value found by HI of 75 ± 15% j3|. 

Note, from Fig. 01 that the perturbative Pomeron contribution to F2 (from scales fx > 
Qo = 1 GeV) is not small; as a rule it is more than half the total contribution. The perturbative 
contribution is even stronger for the 'MRST' fits presented in Table |21 The comparison of the 
separate fits to the ZEUS and HI data presented in Table |21 demonstrates that there is a strong 
correlation between the pairs of parameters q/g and Ci/s, where i = q,g,L. That is, with the 
present accuracy of the data, it is hard to distinguish between partons which originate from 
the two-gluon and two-quark components of the Pomeron.^ Nevertheless, the final diffractive 
parton distributions are similar for the different fits. This stability increases confidence in these 
distributions, so that they can be used in the description of other diffractive processes at HERA 

^In the preliminary HI analysis ^, Aqcd — 0.2 GeV for 4 flavours, corresponding to as{M'^) = 0.1085, 
whereas we take as(Af|) = O.f 190 from the MRST2001 NLO parton set [HI; of. the world average, as(M|) = 
0.1187(20), from the PDG ^. 

^The combined analysis of DDIS data with a more exclusive diffractive process, such as diffractive J/tp 
production at HERA, which is sensitive to the two-gluon component of the Pomeron, may help to resolve this 
problem. 
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Figure 5: (a) The curves labelled 'A' and 'MRST' show the diffractive parton distributions 
extracted from the fits in Tables Q and |21 to the combined preliminary ZEUS ^1 12] and HI 
jH] data (compared to those obtained by HI 0). (b) The breakdown of the five separate 
components of ()20|) for the 'MRST' fit to the combined data sets. 
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and the Tevatron. Of course, we must include the probabihty that the rapidity gap survives 
the soft rescattering of the colhding hadrons or 'hadron-hke' states pT'. 

An advantage of describing the diffractive DIS data using an approach with an exphcit 
integral over the Pomeron scale /i is the possibility to use the results to evaluate the absorptive 
corrections, AFJ''^^(a;B, Q^] fj,'^), to the inclusive structure function F2{xb, Q^)- These corrections 
arise during the DGLAP Ing^ evolution at each point q^ = /i^. With this approach the AGK 
cutting rules [22] give 

AFtix,, Q'; /i^) ^ -F2^(x„ Q'; /i^), (22) 

where F^{xb, Q^\ /i^) is the contribution to the diffractive structure function F^ (integrated 
over xjp) which originates from a perturbative component of the Pomeron of size l//i. In a 
separate paper [23j, we use this property to incorporate the absorptive corrections in a parton 
analysis of inclusive DIS data. 
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